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Abstract. We study the phase diagram of the dilute two-component Fermi gas at zero 
temperature as a function of the polarization and coupling strength. We map out the detailed 
phase separations between superfluid and normal states near the Feshbach resonance. We show 
that there are three different coexistence of superfluid and normal phases corresponding to 
phase separated states between: (I) the partially polarized superfluid and the fully polarized 
normal phases, (II) the unpolarized superfluid and the fully polarized normal phases and (III) 
the unpolarized superfluid and the partially polarized normal phases from strong-coupling BEC 
side to weak-coupling BCS side. For pairing between two species, we found this phase separation 
regime gets wider and moves toward the BEC side for the majority species are heavier but shifts 
to BCS side and becomes narrow if they are lighter. 



1. Introduction: 

Recently experiments in the ultracold Fermi gases [T] have raised strong interest in studying the 
crossover problem from Bose-Einstein condensation (BEC) to the condensation of weak-coupling 
Cooper pairs. The experiments [21 E] on ^Li atoms with imbalance spin populations further 
provide an opportunity to study the superfluid properties with mismatched Fermi surface. The 
phase diagram of this polarized Fermi system has been studied near the Feshbach resonance for 
both of the pairing with equal masses [Ullj and different species [H [3 [H [9l [TO] . The homogeneous 
superfluid is stable only at the large positive coupling strengths. There are two possible scenarios 
to replace the unstable phases, namely phase separation [XT] and FFLO states [T2]. Determining 
phase separation is relatively straight-forward. On the other hand, calculations for the FFLO 
states are highly non-trivial due to the complex spatial dependences of the order parameter 
|13j . However, the instability of the normal state towards the FFLO state is relatively easy 
to determine if one assumes a second order transition, for then one only needs to examine 
the Cooperon at finite q. In this paper, we shall then construct the phase diagram for phase 
separation using a Maxwell construction [14J, but then terminating the procedure whenever the 
corresponding normal state has a finite-q Cooperon instability. The latter regions are presumably 
in the FFLO state. 

For the same species, we found that there are three different coexistence of superfluid 
and normal phases corresponding to the phase separation between (I) the partially polarized 
superfluid (SFp) and the fully polarized normal (Np?) phases, (II) the unpolarized superfluid 
(SFq) and the fully polarized normal phases, and (III) the fully polarized superfluid and the 
partially polarized normal (Npp) states. 



For the pairing between different species, this phase separation region becomes wider 
compared to the pairing of equal masses for the case of the majority species are heavier but 
narrower if the majority species are hghter. For the case with 7 = rrimajority/ ^minority = 6.67 
(between ^'^K and ^Li), the instabihty of the finite-q Cooperon occurs earher than the instabihty 
of the phase separation for almost all couplings so that there are effectively only two kinds of 
phase separation (SFp + Npp and SFo + Npp). For the other case with 7 = 0.15, the phase 
separation region shifts to BCS side and all of three kinds of phase separated states still exist. 



2. Formalism: 

We start from the two-component fermion system across a wide Feshbach resonance which may 
be described by an effective one-channel Hamiltonian as follows: 
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where Co-(k) = h'^k'^ /2m^ — fia-, 9 is the bare coupling strength, and the index a runs over 
the two species {h and /), where —h = I. rrih {mi) to represent the mass of the heavier 
(lighter) component. Within the BCS mean field approximation, the excitation spectrum in 
a homogeneous system for each species is 
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with the reduced mass nir, the chemical potential difference h = (nh — IJ'i)/'^ and the average 
chemical potential /x = {^^ + fJ,i)/2. 

It is convenient to re-scale the energy in unit of the chemical potential difference h, then the 
equation for pairing field A reads 
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where k = k/ y^mr\h\, a = a^Jrr^\h\, and h = 1. The total density (A^) and the polarization (P) 
of the system are 
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with fj, = n/\h\. For a given fl, the scaled scattering length a can be evaluated as a function of 
the pairing field A/|/i| through equation When there are multiple solutions at the same ^ 
and h, the physical solution is then determined by the condition of minimum free energy which 
can be found via the same procedure as the usual Maxwell construction 



3. Results and Discussions 

In Fig. [H we plot the phase diagram of the polarization as a function of the coupling strength 
[{kpa)~^] for the pairing with equal masses. The Fermi momentum (kp = (Sn'^N)^^^) is the 
usual definition for the unpolarized non-interacting Fermi gas. In a polarized system, the 
homogeneous state is stable only at large positive and negative coupling strengths corresponding 
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Figure 1. Phase diagram for the imbalanced 
fermion system with 7 = 1.0. The homogeneous 
states (labehed as SFp and Npp) are stable only 
at large positive or negative coupling strengths. 
The shaded area represents three different phase 
separated states (see text). The short dashed 
lines represents the instability of the finite-g 
Cooperon. The inset shows a typical solution of 
equation ^ at fixed fi in the phase separation 
region. The shaded areas above and below the 
horizontal line are equal. 



to the partially polarized superfluid (BEC side) and the partially polarized normal state (BCS 
side). The shaded area near the resonance (|a| — > 00) represents the phase separation between 
superfluid fluid and normal gas (HE]. Furthermore, this shaded area is divided into three regions: 
the phase separation between (I) SFp and Npp, (II) SFoandNrp, and (III) SFoandNpp. (This 
division is slightly different from the Monte Carlo results [E]). The phase boundaries among 
these three regions are determined by the Maxwell construction at fixed /z and a. A typical 
example is shown in the inset in figure [H The two end points (labelled as N and S) of the 
horizontal line are corresponding to the normal (A = 0) and superfluid (A 7^ 0) states. Then 
the phase boundaries (long dashed lines shown in figured]) are evaluated by 
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where N^'^ = N^'^ + N^'^ , ^d'^ — ^h'^ ~ N^'^ , and < 2; < 1. As increasing fl, we will 
reach the BCS regime, the instabilities of the phase separation and finite q Cooperon intersect 
near P ~ 0.5 and {kpa)"^ ~ —0.49. The FFLO phase probably exists on the right of this point 
L13J. 

For the pairing between different species, the story is quite different. For the case of the 
majority species are heavier, the phase separation region moves toward the BEC side and 
becomes much wider than the case of equal masses pairing. In figure [51 we plot the phase 
diagram for 7 = 6.67. The instability of the finite-g Cooperon occurs earlier compared to the 
instability of the phase separation for almost entire polarization such that the phase separation 
between SFq and Npp effectively does not exist in this case. 

When the majority species are lighter, the phase separation region gets narrow and shifts 
toward the BCS side. In figure [3l we plot the phase diagram for 7 = 0.15. As increasing the 
coupling strengths from BCS limit, the system prefers entering FFLO state first for P < 0.36 
but the phase separated states between the SFoandNpp for P > 0.36. All of the three phase 
separated states exist in this case. 
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4. Conclusion 

We have studied the phase diagram for the pairing between two species in a dilute Fermi gas 
at zero temperature. Near the Feshbach resonance, we map out the detailed phase separations 
between superfluid and normal state. There are three different phase separated states for the 
paring of equal masses. For the majority species are lighter, the phase separation region becomes 
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Figure 2. Phase diagram for the imbalanced 
fermion system with 7 = 6.67. Notations are 
the same as in figure [J but only two kinds of 
phase separated states exist here (see text). 



Figure 3. Phase diagram for the imbalanced 
fermion system with 7 = 0.15. Notations are 
the same as in figure [H 



narrow and shifts to BCS side compared to the case of equal masses. On the other hand, this 
region gets wider and only the phase separated states of SFp + Npp and SFq + Npp exist for 
7 = 6.67. 
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